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METHOD FOR EVALUATING THERMAL COMFORT OF A STRUCTURE AND 
AN ASSISTING METHOD, PROGRAM OR SYSTEM FOR DESIGNING A 
5 STRUCTURE IN CONSIDERATION OF THERMAL COMFORT 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

The present invention relates to a method for 
evaluating thermal comfort of a structure and an 

10 assisting method, program or system for designing a 
structure in consideration of thermal comfort. In 
particular, it relates to a method for evaluating thermal 
comfort of a structure and an assisting method, program 
or system for designing a structure based on 

15 thermophysical properties of a translucent member (such 
as a window glass sheet, a glass plate used as an outer 
wall of an atrium, a plastic plate, a plastic film or the 
like) used for a structure (such as a vehicle, a train, a 
ship, an aircraft, a space shuttle, a building and so on) 

20 and each part (such as a wall or floor material) 

constituting the structure, and based on a shape of each' 
part . 

DISCUSSION OF BACKGROUND 

It has been known that thermal environment in a 
25 vehicle or a building changes largely by the influence of 
solar radiation incident through window glass sheets. In 
order to control such influence of solar radiation, 



various articles such as a solar reflection glass sheet, 
a solar absorption glass sheet, an insulated glazing 
glass sheet and so on are presented by glass 
manufacturers . 

5 These glass sheets are generally evaluated based on 

an insulation performance such as a transmittance of 
solar radiation, a gain of solar heat and so on ruled in 
JIS R 3106. However, when the insulation performance is 
evaluated by only physical characteristics of glass 

10 itself, it is difficult to know thermal sensation that a 
human feels actually. As an attempt of solving such 
problem, there is a proposal that a thermal manikin 
formed in consideration of heat generated from a human 
body (a dummy doll for reproducing thermal 

15 characteristics of a human) is prepared and a heat loss 
from each human part (such as a head part, a chest part, 
arms, legs or the like) is measured by using the thermal 
manikin to examine the thermal sensation. 

However, even when such thermal manikin is used, it 

20 is difficult to examine individually a thermal effect 

that each part (a window glass sheet, a seat or the like) 
constituting a vehicle body affects a passenger. In the 
conventional technique, there is no way other than the 
study of the influence obtained by considering all 

25 factors of solar irradiation transmitting through window 
glass sheets, and radiation and convection on each part. 
Accordingly, it was difficult to examine quantitatively 
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contribution of each part, and it was difficult to select 
the material, property or shape of each part which can 
provide the optimum thermal comfort in the passenger 
compartment . 

5 Further, even though the thermal effect of each part 

could be examined by using the thermal manikin, an actual 
examination was not easy because a prototype of a • room or 
the passenger compartment of a vehicle was required in 
order to evaluate the thermal environment of the room or 

10 passenger compartment. Further, the thermal manikin is a 
very expensive device because it is provided with a large 
number of sensors and heat sources to reproduce the 
thermo-regulating function of a human. A patent document 
1 ( JP-A-2001-155052 , p. 1-7) discloses a method for 

15 selecting a glass sheet and a patent document 2 (JP-A-10- 
239161, p. 1-8) discloses a technique of simulating 
thermal environment. However, these documents do not 
disclose that simulation is conducted in consideration of 
a precise shape of human body. 

20 It is an object of the present invention to provide 

a method for evaluating thermal comfort of a structure 
and an assisting method, program or system for designing 
a structure in consideration of thermal comfort, which 
can easily be carried out without requiring a laboratory 

25 equipment such as a prototype. 

SUMMARY OF THE INVENTION 
The present invention is provided in order to 



achieve the above-mentioned object. 

In accordance with a first aspect of the present 
invention, there is provided a method for evaluating 
thermal comfort of a structure constituted by a first, 
5 second, ... n-th part (n is a natural number of at least 

two) wherein at least one of these parts is a translucent 
member for introducing light to the inside, the method 
for evaluating thermal comfort of a structure comprising: 
a step (a) of preparing data of (1) the material 

10 property of the first part, the material property of the 
second part, ... the material property of the n-th part, 
and data of (2) an mount of solar radiation passing 
through the translucent member to reach a measuring 
device having a shape imitating a human body part, an 

15 amount of solar radiation to the structure, an amount of 
convection heat transfer in the structure, an amount of 
radiation heat transfer in the structure, humidity in the 
structure and/or a thermo-regulating function of the 
measuring device, and calculating at least one of the 

20 amount of heat loss from the surface of the measuring 
device, the temperature of the measuring device and/or 
the wettedness at the surface of the measuring device 
based on at least one in each of data (1) and (2), and 

a step (b) of calculating a thermal comfort index of 

25 the measuring device by using a result of the above 
calculation . 

According to a second aspect of the present 
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invention, there is provided an assisting method for 
designing a structure in consideration of thermal 
comfort, the structure being constituted by a first, a 
second, ... an n-th part (n is a natural number of at least 
5 two) wherein at least one of these parts is a translucent 
member for introducing light to the inside, the assisting 
method for designing a structure comprising: 

a step (a) of preparing at least one candidate value 
for each of the material property of the first part, the 

10 material property of the second part, ... the material 
property of the n-th part, and registering candidate 
values of material properties of these parts in a 
predetermined data base, 

a step (b) of selecting and reading the material 

15 property registered in said data base for each part, and 
preparing at least one combination among the registered 
material properties, 

calculating at least one of the amount of heat loss 
from the surface of the measuring device, the temperature 

20 of the measuring device and/or the wettedness at the 

surface of the measuring device based on at least one in 
each of data (1) combinations of material properties thus 
prepared and (2) an mount of solar radiation passing 
through the translucent member to reach the measuring 

2 5 device having a shape imitating a human body part, an 

amount of solar radiation to the structure, an amount of 
convection heat transfer in the structure, an amount of 
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radiation heat transfer in the structure, humidity in the 
structure and/or a thermo-regulating function of the 
measuring device, and 

calculating a thermal comfort index of the measuring 
5 device by using a result of the above calculation, and 
a step (c) of selecting a combination of material 
properties having a value closest to the optimum thermal 
comfort index . 

According to a third aspect of the present 

10 invention, there is provided a thermal comfort evaluation 
program for a structure constituted by a first, second, ... 
n-th part (n is a natural number of at least two) wherein 
at least one of these parts is a translucent member for 
introducing light to the inside, the thermal comfort 

15 evaluation program comprising program codes for 

instructing to a computer to realize the steps described 
in the first aspect. 

According to a fourth aspect of the present 
invention, there is provided a thermal comfort evaluation 

20 system for a structure constituted by a first, a second, 
... an n-th part (n is a natural number of at least two) 
wherein at least one of these parts is a translucent 
member for introducing light to the inside, the thermal 
comfort evaluation system for a structure comprising: 

25 means for preparing data of (1) the material property 

of the first part, the material property of the second 
part, ... the material property of the n-th part, and data 



of (2) an mount of solar radiation passing through the 
translucent member to reach a measuring device having a 
shape imitating a human body part, an amount of solar 
radiation to the structure, an amount of convection heat 
transfer in the structure, an amount of radiation heat 
transfer in the structure, humidity in the structure 
and/or a thermo-regulating function of the measuring 
device, and calculating at least one of the amount of 
heat loss from the surface of the measuring device, the 
temperature of the measuring device and/ or the wettedness 
at the surface of the measuring, device based on at least 
one in each of data (1) and (2), and 

means for calculating a thermal comfort index of the 
measuring device by using a result of the above 
calculation. 

According to a fifth aspect of the present 
invention, there is provided an assisting program for 
designing a structure constituted by a first, second, ... 
n-th part (n is a natural number of at least two) wherein 
at least one of these parts is a translucent member for 
introducing light to the inside, the assisting program 
for designing a structure comprising program codes for 
instructing to a computer to realize the steps described 
in the second aspect. 

According to a sixth aspect of the present 
invention, there is provided an assisting system for 
designing a structure constituted by a first, a second, ... 



an n-th part (n is a natural number of at least two) 
wherein at least one of these parts is a translucent 
member for introducing light to the inside, the assisting 
system for designing a structure comprising: 
5 (a) means for preparing at least one candidate value 

for each of the material property of the first part, the 
material property of the second part, ... the material 
property of the n-th part, and registering candidate 
values of material properties of these parts in a 

10 predetermined data base, 

(b) means for selecting and reading the material 
property registered in said data base for each part and 
preparing at least one combination among the registered 
material properties , 

15 calculating at least one of the amount of heat loss 

from the surface of the measuring device, the temperature 
of the measuring device and/or the wettedness at the 
surface of the measuring device based on at least one in 
each of data (1) combinations of material properties thus 

20 prepared and (2) an mount of solar radiation passing 
through the translucent member to reach the measuring 
device in a form of a human body part, an amount of solar 
radiation to the structure, an amount of convection heat 
transfer in the structure, an amount of radiation heat 

25 transfer in the structure, humidity in the structure 
and/or a thermo-regulating function of the measuring 
device, and 



calculating a thermal comfort index of the measuring 
device by using a result of the above calculation, and 

(c) means for selecting a combination of material 
properties having a value closest to the optimum thermal 
comfort index. 

BRIEF DESCRIPTION OF THE DRAWINGS 
A more complete appreciation of the invention and 
many of the attendant advantages thereof will be readily 
obtained as the same becomes better understood by 
reference to the following detailed description when 
considered in connection with the accompanied drawings, 
wherein: 

Fig. 1 is a block diagram showing an embodiment of 
the thermal comfort evaluation apparatus according to the 
present invention ; 

Fig. 2 is a partly broken perspective view showing a 
numerical model to evaluate the thermal comfort of a 
window glass sheet; 

Fig. 3(a) is a partly broken perspective view 
showing each part constituting a vehicle body (viewed 
from a front side) and Fig. 3(b) is a partly broken 
perspective view showing each part constituting a vehicle 
body (viewed from a back side) ; 

Fig. 4(a) is a partly broken perspective view 
showing the behavior of direct solar radiation in Type I, 
and Fig. 4(b) is a similar view in Type II; 

Fig. 5 is a diagram showing a thermal model of human 



body shape; 

Fig. 6(a) is a front view showing classification of 
each part of a surface model of human body, Fig. 6(b) is 
a front view showing the detail (the arrangement of 
surface elements) of the surface model of human body, and 
Fig. 6(c) is a side view showing the detail (the 
arrangement of surface elements) of the surface model of 
human body; 

Fig. 7 is a block diagram showing a thermo- 
regulation model; 

Fig. 8 is a flow chart showing a sequence of 
selection (main routine) of a window glass sheet and each 
part; 

Fig. 9 is a flow chart showing the detail of Step 

104; 

Fig. 10 is a flow chart showing the detail of Step 

201; 

Fig. 11 is a diagram showing various kinds of solar 
radiation ; 

Fig. 12 is a graph showing the relation between an 
incident angle and transmittance as well as reflectivity; 
Fig. 13 is a flow chart showing the detail of Step 

309; 

Fig. 14 is a partly broken perspective view showing 
an analytical space and a surface model of thermal 
manikin in a case that a vehicle of recreational type 
accommodates a passenger; 



Fig. 15(a) is a partly broken perspective view 
showing a result of simulation using solar radiation and 
Fig. 15(b) is a partly broken perspective view showing a 
result of simulation using a climate control room; 
5 Fig. 16 is a graph showing an amount of absorptance 

of solar radiation received by the entirety of a manikin; 

Fig. 17 is a graph showing a heat load by an air 
conditioner; and 

Fig. 18 is a graph showing a result of numerical 
10 experiments. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the following, preferred embodiments of the 
present invention will be described with reference to 
drawings . 

15 Fig. 1 is a block diagram showing an embodiment of 

the present invention. In Fig. 1, a thermal comfort 
evaluation apparatus 10 comprises a computer system such 
as a work station or the like, which is connected to 
external devices such as HDD (Hard Disk Drive) , an 

20 external memory unit 21 such as an optical disk device, 
an electro-optical disk device or the like, a 
manipulating part 22 such as a key board, a mouse or the 
like, and a display device 23 such as CRT (Cathode Ray 
Tube) , a liquid crystal display, PDP (Plasma Display 

25 Panel) or the like. 

The thermal comfort evaluation apparatus 10 
comprises a central processing unit 11, a memory part 12 
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such as RAM (Random Access Memory) , or ROM (Read Only 
Memory) or the like, I/F13, 14 and 15 functioning as 
input/output interfaces and a bus 16. The bus 16 is 
connected to each part in the thermal comfort evaluation 
5 apparatus 10 so that each part sends or receives an 
address signal, a data signal and various kinds of 
control signals through the bus 16. 

The central processing unit 11 is a device having an 
operating function and a controlling function, and it 

10 reads out a program for a numerical simulation from the 
external memory unit 21 to execute operation in the 
memory part 12 as a working area. Accordingly, the 
central processing unit 11 is provided with an operation 
part 11a, an evaluation part lib and a selection part 11c 

15 by means of hardware or software to execute a numerical 
simulation . 

The manipulation part 22 is connected to an I/F13 . 
When a user operates the manipulation part 22, an 
instruction signal is supplied through the I/F13 to the 

20 central processing unit 11. The external memory device 

21 which holds a database including various kinds of data 
such as data of the shape of structures and human bodies, 
data concerning calculation conditions, data of materials 
and shapes (which are specified by material properties) 

25 usable for each part of a vehicle body (such as main 
body, seat, instrument panel, rear panel, toe board, 
center console and so on) and many types of glass, is 



- 13 - 

connected to an I/F14. The reading-out of data from the 
external memory device 21 and the writing-in of data to 
the external memory device 21 are carried out under 
control by the central processing unit 11. I/F15 is 
5 connected to the display device 23. In the display 

device 23, information inputted by the manipulating part 
22 or a result simulation is visually displayed. 

The thermal comfort evaluation apparatus 10 reads 
out data of shapes of a structure or a numerical model of 

10 human body from the external memory device 21 to perform 
a simulation of thermal condition in an indoor space or a 
passenger compartment by using the obtained data. A 
result of the simulation is stored in the external memory 
device 21 and a visual image is displayed in the display 

15 device 23. Accordingly, a user can evaluate the comfort 
of a translucent member such as a window glass sheet 
attached to a structure or a vehicle by observing such 
image. For instance, the evaluation of thermal 
environment in a passenger compartment is carried out as 

20 follows. 

Fig. 2 is a partly broken perspective view showing 
an example of evaluating the thermal environment of a 
passenger compartment. Fig. 2 shows a surface model 
comprising a plurality of surface elements of a passenger 
25 compartment (a sedan-type vehicle of right starring wheel 
accommodating 5 passengers) and passengers in the 
compartment. Two separated seats are provided in a front 
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portion of the passenger compartment and a bench type 
seat capable of accommodating 3 passengers is provided in 
a rear portion of it. In the surface model, 4 passengers 
in total are on the front and rear seats and the 
5 passenger on the driver's seat holds the starring wheel. 

A windshield is provided at a front portion of the 
compartment, sidelites (a front sidelite and a rear 
sidelite) are provided at both sides and a rear glass is 
provided at a rear portion of the compartment. Solar 

10 radiation over a window glass sheet attenuates in 

response to material properties of glass (i.e., the 
transmittance of solar radiation, the absorptance of 
solar radiation, the emissivity and the overall heat 
transfer coefficient) before entering into the passenger 

15 compartment. The incident light increases compartment 
temperature because it heats each part in the passenger 
compartment . 

On the other hand, surfaces of the passengers and 
each vehicle interior part are divided into a plurality 

2 0 of surface elements , and the indoor space between the 

passengers and vehicle interior parts is divided into a 
plurality of cubic elements (not shown) . In this 
embodiment, a quadrangle surface element expressed by a 
bilinear form and a hexahedra element are used. By using 

25 these surface elements and cubic elements, a numerical 
simulation about the solar radiation to the passenger 
compartment , radiation in the compartment and temperature 
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and airflow fields in the indoor space is carried out. 

In the passenger compartment, it is unavoidable that 
the solar radiation entering through window glass sheets 
provides a thermally non-uniform environment. On the 
5 assumption that a heating, ventilation and air 

conditioning system (not shown) is installed in the 
passenger compartment, positions of an airflow outlet and 
an airflow inlet, temperature of flowing air, rates of 
flowing air and sucked air, directions of flowing air and 

10 sucked air, the absolute humidity (or relative humidity) 
of air and so on can be determined optionally. Data 
relating to this heating, ventilation and air 
conditioning system are registered in the external memory 
device 21 in Fig. 1 as well as material properties of the 

15 parts. 

Other than the heating, ventilation and air 
conditioning system, a radiation cooling/heating panel 
system, a ventilation or a humidity control system may be 
provided optionally. Calculation conditions used in this 

20 simulation are shown in Table 1. Further, since the 

position of the sun can theoretically calculated from the 
present location of the vehicle and date and hour on 
which the simulation is carried out, it is possible to 
determine an irradiation angle of solar radiation 

2 5 entering into the passenger compartment through the 
window glass sheets. 
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Table 1 Calculation conditions 

Present location : Tokyo 

Date and hour: 12:00, August 1 

Outdoor temperature: 3 5°C 

Human body: (metabolic rate: 1 met, clothing amount: 
0.6 clo, relative humidity: 50%, 
absorptance of solar radiation: 0.74) 
Temperature of flowing air at a front seat : 10°C 
Running forward the south at 40 km/h 

Fig. 3 is a partly broken perspective view showing 
each part constituting a vehicle body. Fig. 3(a) is such 
5 view from a viewpoint looking down the windshield from a 
front side of the vehicle , and there are provided the 
windshield and sidelites (front and rear) in the vehicle 
body. A ceiling is provided from an upper edge of the 
windshield to a rear portion of the vehicle, and a rear 

10 glass sheet is disposed in the rear portion of the 
vehicle. Fig. 3(b) is a view looking the passenger 
compartment from a rear portion of the vehicle, and parts 
such as front and rear seats, a rear panel, a toe board, 
an instrument panel, a front panel, a center consol and 

15 so on are provided on the shasses (not shown) . 

Fig. 4 is a diagram showing a distribution of direct 
solar radiation incident into a passenger compartment in 
a case that two types (Type I and Type II) shown in Table 
2) are used. Fig. 4(a) shows the case that glass of Type 

20 I is used, and Fig. 4(b) shows the case that glass of 
Type II is used. The intensity in a distribution of 
solar radiation is expressed by gradation of white and 
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black in Fig. 4 (namely, the solar radiation receiving 
quantity is larger as a white tone increases, and the 
solar radiation receiving quantity is smaller as a black 
tone increases) . As understood from the observation to 
the passenger on the driver seat, the amount of direct 
solar radiation reaching the passenger compartment is 
smaller in Type II and therefore, Type II is excellent in 
insulation performance. 



Table 2 (a) Type I 





Type of glass 


Transmit tance 
of solar 
radiation (%) 


Absorp tance 
of solar 
radiation (%) 


Windshield 


Glass A 


52 .7 


41.7 


Sidelite 
(front) 


Glass B 


60.0 


34 . 0 


Sidelite 
(rear) 


Glass B 


60 . 0 


34 . 0 


Rear glass 


Glass B 


60 . 0 


34.0 



(b) Type II 





Type of glass 


Transmit tance 
of solar 
radiation (% ) 


Absorp tance 
of solar 
radiation (%) 


Windshield 


Glass C 


45.3 


28 . 7 


Sidelite 
(front) 


Glass D 


48.4 


46.1 


Sidelite 
(rear) 


Glass E 


18.4 


77 . 1 


Rear glass 


Glass E 


18.4 


77 . 1 



Further, the passenger compartment is constituted by 
n (n is a natural number of at least two) number of parts 
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including at least a translucent member, and the material 
property of the material used for each part is stored in 
the external memory device 21 in Fig. 1 as well as 
properties of two types of glass in Table 2 . The 
5 material property usable for each part can be shown in a 
matrix form as shown in Table 3(a). For example, with 
respect to a part P, any one of material properties M 0 o , 
Mio, M 2 o/ M 30/ M 4 o and M 50 can be selected. 

Namely, at least one material property M** (=M 0 o# M 0 i, 

10 Mio, ...) is registered for each part, and each material 

property correspond to a candidate value of a material 
for an actually used part. When an instrument panel is 
taken as an example , material properties of plastic 
materials are chosen and when a seat is taken as an 

15 example, material properties of fabrics or leather are 
chosen. Further, the material property M** includes at 
least one selected from the group consisting of 
transmittance of solar radiation, absorptance of solar 
radiation, emissivity , overall heat transfer coefficient 

20 and data (CAD data) of each part. 

Accordingly, it is necessary to prepare at least one 
combination of material properties M** of parts 
constituting the passenger compartment after the material 
property M** is selected for each part (Table 3 (b) ) . 

25 Then, a thermal comfort index that a human feels (SET*, 
equivalent temperature or TSV or the like, which will be 
described later) is calculated for each combination. By 
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comparing thermal comfort indices obtained by 
calculation, it is possible to obtain a combination of 
material properties which can provide the optimum thermal 
comfort as shown in Table 3(c). Here, a value closest to 
5 "0" is an optimum value. 

In this case, when the thermal comfort index is 
calculated in consideration of not only the material 
property but also data relating to a heating, ventilation 
and air conditioning system or the like, it is possible 

10 to design an optimum structure with such heating, 

ventilation and air conditioning system. In this case 
also, data relating to the air-conditioning system are 
read out from the external memory device 21 in Fig. 1 as 
well as the data of material properties, and the thermal 

15 comfort index is calculated by adding the read-out date. 
Further, when data concerning various shapes of each part 
are prepared as in Table 3, it is possible to select a 
combination of shapes of parts, which provides an optimum 
thermal comfort . 



20 
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Table 3(a) Material property usable for each part 







Part P 


Part Q 


Part R 


Part S 


Material 
property 


M 00 


M01 


M02 


M 03 


M 10 


Mn 


M12 


M 13 


M 2 o 


M21 


M 22 


M 23 


M 30 


M31 


M 32 


M33 ! 




M 40 


M41 


M 42 


M43 


M50 


M51 


M 52 


M53 


M** " ~ — — 




Transmittance of solar radiation, 
absorptance of solar radiation, 
emissivity, overall heat transfer 
coefficient etc. 





Table 3 (b) Example of combination of material properties 



Part P 


Part Q 


Part R 


Part S 


TSV 


M10 


M 2 i 


M12 


M 23 


2 (warm) 


M10 


M31 


M 32 


M 33 


3 (hot) 




M30 


M31 


M 52 


M43 


0 (neutral) 





5 Table 3(c) Example of combination capable of obtaining 



the optimum thermal comfort index 



Part P 


Part Q 


Part R 


Part S 


M10 


M31 


M 52 


M 43 



Next, description will be made as to a thermal model 
of human body in consideration of solar radiation, 
10 proposed by the present invention. 

Fig. 5 is a diagram showing a thermal model of human 
body. A manikin for evaluating thermal environment 
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prepared based on the thermal model is a so-called 
thermal manikin. The thermal manikin used in the present 
invention is called a numerical thermal manikin which can 
reproduce thermal sensation of human body by only 
5 computer simulation. 

There have been proposed several thermo-regulation 
models of human body. As representative examples of 
these proposals, there are a two-node model by A. P. 
Gagge, A. P. Fobelets' and L.G. Berglund (see a Standard 

10 Predictive Index of Human Response to the Thermal 

Environment, ASHRAE Transactions, Vol. 92, pp. 7 09-731, 
1986.) and a 65 multi-nodes model by S. Tanabe, et al 
(see EVALUATION OF THERMAL COMFORT USING COMBINED MULTI- 
NODE THERMOREGULATION (6 5MN) AND RADIATION MODELS AND 

15 COMPUTATIONAL FLUID DYNAMICS (CFD) , Energy & Buildings- 
Special Issue, 2002. 4). 

However, the above-mentioned proposals do not take 
the shape of human body sufficiently into consideration. 
For example, in the two-nodes model, a model of human 

20 body is formed by a core layer and a shell layer covering 
the core layer. Accordingly, simulation is impossible 
with respect to solar radiation to a part such as the 
head, a hand or a foot. The 65 multi-nodes model has a 
more precise structure of human body than the two-nodes 

25 model, namely, it has 16 parts such as the head, chest, 
legs and so on and the surface area and the weight of 
each part are specified. However, consideration is not 
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taken sufficiently with respect to the shape of each 
part. Since the shape of each part is unclear even in 
the 65 multi-nodes model, the angle factor between each 
part and a wall or a window can not be calculated, and 
5 therefore, it is difficult to obtain correctly the 
distribution of solar radiation to the human body. 

On the other hand, there have been proposed several 
simulation systems in consideration of the shape of human 
body. However, a sufficient study has not been made on 

10 environmental conditions around a human body (solar 
radiation, heat radiation or the like) in order to 
simulate thermal sensation of human body correctly (see 
G. Fronzes, et al . : Athebes, a car thermal behavior 
simulation software, ATA-Ingegneria Automotoristica, 

15 Marzo, Vol. 49, No. 3, 1996.). 

Accordingly, it is difficult to simulate correctly 
thermal sensation in conditions that the passenger 
compartment is influenced largely by solar radiation. 
According to the present invention, the above- 

2 0 mentioned problem can be eliminated by combining a model 
having a shape imitating the shape of a human body 
(hereinbelow, referred to as a surface model of human 
body or simply, a surface model) and a model having a 
human thermo-regulation function (hereinbelow, referred 

25 to as a thermo-regulation model) , whereby the thermo- 
regulating function of a human body can correctly be 
simulated even in an environment that influence by solar 
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radiation is large. 

From the above-mentioned viewpoint, a thermal model 
of human body (a thermal manikin) 3 0 shown in Fig. 5 is 
formed by the combination of a surface model 31 and a 
5 thermo-regulation model 32. Figs. 6 and 7 show 

particulars of the structure. Fig. 6(a) is a front view 
showing roughly classified parts of the surface model 
shown in Fig. 5, and Figs. 6(b) and 6(c) are respectively 
a front view and a side view showing the detail (surface 

10 elements) of the surface model. 

The model shown in Fig. 6(a) is classified into 16 
parts in total: head 31-1, chest 31-2, back 31-3, pelvis 
31-4, 1-shoulder 31-5, r-shoulder 31-6, 1-arm 31-7, r-arm 
31-8, 1-hand 31-9, r-hand 31-10, 1-thigh 31-11, r-thigh 

15 31-12, 1-leg 31-13, r-leg 31-14, 1-foot 31-15, and r-foot 
31-16. However, human parts is not necessary to be 
classify into 16 parts but may be changed optionally. 

In order to simplify simulation of thermo-regulating 
function of each part, it can be also assumed that each 

20 part has a skin layer installing therein a heat source (a 
clothing layer is provided depending on a part) , and the 
thermo-regulation model shown in Fig. 7 is provided for 
each part. The surface area and the weight of each part 
are shown in Table 4 . 



25 
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Table 4 Parts of 65 multi-nodes model 



Body part 


Surface area of 
body part (m 2 ) 


Weight of body 
part (kg) 


1 . Head 


0 . 140 


4 . 020 


2 . Chest 


0 . 175 


12 .400 


3 . Back 


0 . 161 


11.030 


4. Pelvis 


0 . 221 


17 . 570 


5. L- shoulder 


0 . 096 


2 . 163 


6. R- shoulder 


0.096 


2 . 163 


7 . L- arm 


0 . 063 


1.373 


8. R-arm 


0 .063 


1 .373 


9 . L-hand 


0 .050 


0.335 


10. R-hand 


0.050 


0 . 335 


11. L-thigh 


0.209 


7 . 013 


12. R- thigh 


0 .209 


7 . 013 


13. L-leg 


0 .112 


3 .343 


14. R-leg 


0 . 112 


3 .343 


15. L-foot 


0 .056 


0 . 480 


16. R-foot 


0 .056 


0 .480 


The whole body- 


1.870 


74.430 



As shown in Figs. 6(b) and 6(c), the surface of each 
part is divided into a plurality of surface elements, and 
5 calculation of an amount of solar radiation reaching the 
human body, an amount of heat loss radiated from the 
human body surface, skin temperature and so on is carried 
out for each surface element. 

Next, description will be made as to a thermo- 
10 regulation model assembled in each part of a human 
figure . 

Fig. 7 is a block diagram showing the thermo- 
regulation model assembled in each part. A thermo- 
regulation model 32 comprises 16 parts (the k-th part 34 
15 comprises a skin layer 34d including therein a heat 
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source 34d-l and a temperature sensor 34d-2 and a 
clothing layer 34e) . A heat quantity generated from the 
heat source 34d-l is controlled by a personal computer 
3 4d~3. A heat transport between the skin layer and the 
5 clothing layer is effected as follows. In the skin layer 
34d, the heat quantity from the heat source 34d-l is 
controlled based on temperature measured by the 
temperature sensor 34d-2 so that the skin temperature is 
controlled. There is a heat transport due to thermal 
10 conduction between the clothing layer 34e and the skin 
layer 34d. There is a case that an air layer exists 
between the clothing layer and the skin layer. 
Accordingly, in conducting a simulation, it is necessary 
to consider a heat transport due to any one among 
15 convection, radiation and humidity. The heat generation 
means may be disposed in a layer other than the skin 
layer (e.g., inside the human figure). 

Since the clothing layer 34e is exposed to an 
external environment 40, there is a heat transport due to 
20 convection and radiation. On the assumption that there 
is a heat source 41 such as the sun or the like (or, a 
lamp may be disposed artificially) in the external 
environment 40, the clothing layer 34e is warmed by 
radiation from the heat source 41 (or light beams 
25 irradiated by the lamp) . The presence or absence of the 
clothing layer 34e can be determined for each surface 
element. In this embodiment > calculation is conducted 
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under conditions that the skin layer is exposed at the 
head 31-1, the 1-hand 31-9 and the r-hand 31-10, and 
therefore they have no clothing layer. Further, the 
wavelength of light beams irradiated from the heat source 
5 41 can be optionally determined by the manipulating part 
22, and set values are stored in the external memory 
device 21 in Fig. 1. 

Next, description will be made as to a sequence of 
evaluating the thermal comfort of a passenger caused by a 

10 window glass by a numerical simulation, and selecting the 
optimum type of glass as a window glass based on a result 
of the evaluation, and selecting a combination of parts 
capable of providing the optimum thermal environment. 
Fig. 8 is a flow chart showing a sequence of 

15 selecting material properties of glass and each part. 
First, CAD (Computer Aid Design) data indicating a 
previously prepared vehicle body shape, a plurality of 
surface elements constituting an inner surface of the 
passenger compartment, CAD data indicating a human body 

20 figure, a plurality of surface elements indicating shapes 
of the human body surface, a plurality of cubic elements 
obtained by dividing the indoor space between parts 
constituting the passenger compartment and human body 
figure, and candidate values of some types of glass 

25 usable for window glass sheets (in this embodiment, Types 
I and II shown in Table 2) are stored in the external 
memory device 21 in Fig. 1 (Steps 101 and 102). Further, 
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the material property of each part as shown in Table 3 is 
stored in the external memory device 21 (Step 103). 

Then, the operation part 11a and the evaluation part 
lib read out data of types of glass and their material 
5 properties stored in the external memory device 21 to 
prepare possible combinations of these data, and 
calculation of solar radiation passing through window 
glass sheets, calculation of the thermo-regulation model 
of human body figure and calculation of heat transport in 

10 the passenger compartment (i.e. calculations of 

convection, radiation, humidity and so on) are conducted 
in a combined manner to thereby perform a numerical 
simulation of thermal sensation that a human is expected 
to feel (Step 104) . Alternatively, a simulation may be 

15 carried out by changing material properties in response 
to angle factors shown in Fig. 14 as described later. 

Then, determination is made as to whether or not 
thermal comfort indices (equivalent temperature, a 
predicted mean vote (PMV) or a thermal sensation vote 

20 (TSV) ) are calculated with respect to all combinations 
every time the simulation is completed (Step 105) . If 
not, the operation is returned to Step 104. After the 
calculation has been made with respect to all 
combinations of types of glass and material properties of 

25 parts, the selection part 11c in Fig. 1 compares each 

result of simulation, and selects a combination of a type 
of glass and the material property of each part which 
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brings the best result (Step 106) . As a result, the 
combination of material properties as shown in Table 3(c) 
can be selected for example. When it is enough to 
evaluate only the thermal comfort of a vehicle body 
5 (structure) , it is sufficient to conduct operations from 
the beginning to Step 105. 

There are some options in selecting the combination 
which brings the best result. There can be considered to 
select "neutral" in terms of TSV, which will be described 

10 thereafter, in an intermediate season, to select 

"slightly cool" in summer and to select "slightly warm" 
in winter. However, the determination of these options 
depends on a performer of the simulation. As a possible 
method, there is considered a method that the optimum 

15 value is previously determined, and a value closest to 
the optimum value is selected. 

Further, in selecting glass having excellent thermal 
performance and excellent translucent performance, 
several kinds of glass having excellent thermal comfort 

20 index are selected, and then, glass having the best 
visible light permeability should be selected. 
Combinations of several types of glass usable for window 
glass sheets are not limited to those shown in Table 2 
but it is possible to determine combinations optionally. 

25 Description will be made in detail with reference to 

Figs. 9 to 13 concerning a numerical simulation of 
thermal sensation . 
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Fig. 9 is a flow chart showing particulars of Step 
104. In this embodiment, first, an equivalent 
temperature as an index indicating thermal sensation is 
calculated (Step 201) , and TVS is calculated by using the 
5 equivalent temperature (Step 202). TSV is an index 
corresponding to the thermal sensation that a human 
actually feels, which is calculated according to a 
conversion formula in Table 5 in this embodiment. 

There is a difference in thermal sensation depending 

10 on seasons and areas. Table 5 shows regression equations 
in Japan (summer season), Japan (autumn season), U.S.A., 
Denmark and Singapore. Table 6 shows correspondence 
between thermal sensation votes (TSV) and thermal 
sensations that a human feels. Table 6 shows 7 ranks 

15 wherein "0" means "neutral", "1" slightly warm, "2" warm, 
"3" hot, slightly cool, "-2" cool and "-3" cold. 

Temperature that thermal sensation vote of human is 
neutral (neutral temperature) is described adjacent to 
regression equations . 

20 

Table 5 





Regression equation 


Neutral 
temperature ( °C ) 


Japan 
(summer season) 


TSV=-8 .8 82+0 . 339 -SET* 


26.2 


Japan 
( autumn season ) 


TSV=-6 . 883+0.271- SET* 


25.4 


U.S.A. 


TSV=- 8. 010+0 .290 -SET* 


25.8 


Denmark 


TSV=-7 . 42 0+0 . 310 • SET* 


25 . 6 


Singapore 


TSV=-9 . 3 88+0 . 3 67 -SET* 


25.6 



Table 6 



-3 



-2 



-1 



0 



1 



2 



3 




► 



Cold 



Cool 



Slightly 
cool 



Neutral 



Slightly 



Warm 



hot 



warm 



Fig. 10 is a flow chart showing particulars of Step 
201. First, as a preparatory operation of calculation, 
operations from Step 301 to Step 304 are carried out. 
Then, the inner surface of the passenger compartment and 
the shape of surface of thermal manikin are divided into 
a plurality of surface elements and the indoor space 
between the inner surface of the passenger compartment 
and the thermal manikin are divided into a plurality of 
cubic elements (Step 301) . Then, calculation of angle 
factor is made with respect to all surface elements of 
the thermal manikin and the inner surface of the 
passenger compartment (Step 3 02) . The angle factor is a 
parameter of non-dimensional number which can determine 
radiation exchange between surface elements. 

Then, thermal conditions such as the overall heat 
transfer coefficient, the emissivity, the absorptance of 
solar radiation, the transmittance of solar radiation and 
so on are assigned to each surface element (Step 3 03) . 
Calculation conditions, shown in Table 1, such as an 
amount of thermal insulation of clothes of thermal 
manikin, a present location of the vehicle (i.e., a 
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latitude and a longitude) and date and hour and so on are 
determined (Step 3 04) . These calculation conditions are 
inputted by means of the manipulating part 22 in Fig. 1, 
and inputted data are stored in the external memory 
5 device 21. 

Then, an azimuth and an altitude at the position of 
the vehicle are calculated based on the said latitude and 
a longitude and the season in the operation part 11a so 
that the irradiation angle of solar radiation received by 

10 the vehicle is calculated. Then, various kinds of solar 
radiation quantity reaching directly a human body (any 
one of quantities of direct solar radiation, sky-diffused 
solar radiation, solar radiation reflected on the ground 
and internal reflection of the solar radiation or an 

15 optional combination of these) are calculated based on 
the irradiation angle and properties of window glass 
sheets (Table 2) (Step 305) . Then, the amount of 
absorptance of solar radiation is calculated based on the 
above-mentioned result of calculation and the absorptance 

20 of clothing surface or skin surface. 

As described in the document 1, solar radiation 
includes light beams of a wavelength of 0.3-2.5 \im. 
Accordingly, it is necessary to consider the wavelength 
in calculating the solar radiation quantity. For 

2 5 example, calculation may be conducted over the entire 
wavelength region of 0.3-2.5 urn or calculation may be 
conducted for each wavelength by selecting optionally 
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some kinds of wavelengths. Then, calculation of heat 
transfer in the structure and calculation of the thermo- 
regulation reaction are conducted in combination (Steps 
306-308) and the equivalent temperature at the skin 
5 surface of the thermal manikin is calculated (Step 3 09) . 
Description will be made as to particulars of 
calculation of solar radiation. An indoor space defined 
by a wall member containing a translucent member 
undergoes largely influence of transmitting solar 

10 radiation. Accordingly, accuracy in the prediction of 

thermal comfort depends largely on how a distribution of 
solar radiation gain from a heat source can be predicted 
correctly. It is unnecessary to conduct the calculation 
of solar radiation together with calculation of 

15 temperature and airflow but it can be conducted 

independently. For example, after calculations of an 
amount of direct solar radiation reaching a surface of 
thermal manikin and a surface of the inner wall of the 
indoor space, an amount of sky-diffusion solar radiation, 

20 an amount of solar radiation reflected on the ground and 
an amount of relative reflection of solar radiation 
inside the reached solar radiation, the amount of solar 
radiation gain at the surface of the thermal manikin and 
the surface of the inner walls is calculated (Fig. 11) on 

25 the basis of each solar radiation mentioned above. Then, 
temperature and airflow are calculated on the basis of 
the amount of solar radiation gain at the surface of the 



- 33 - 

thermal manikin and the surface of the inner wall. 
Further, the temperature inside the thermal manikins may 
be calculated. 

The sequence of calculations is as follows. First, 
5 the position of the sun is calculated by inputting the 

present location of the building and date and hour at the 
present time. Then, an amount of direct solar radiation 
Idn to the normal plane and an amount of sky-diffusion 
solar radiation I sky to the horizontal plane are estimated 

10 by an empirical equation based on the obtained data of 
measurements. Then, various amounts of reached solar 
radiation are calculated based on the thus obtained 
amount of direct solar radiation to the normal plane and 
the amount of sky-diffusion solar radiation to the 

15 horizontal plane by taking the thermal performance of the 
wall such as the transmittance of solar radiation t± of a 
translucent member, the reflectivity pi and so on, and 
geometrical shapes of parts in the indoor space and a 
solar radiation shading member such as a hood, into 

20 consideration . 

A method for calculating an amount of solar 
radiation reaching the surface of a thermal manikin is as 
follows . 

1) Calculation of an amount of direct solar 
25 radiation 

An amount of direct solar radiation I di reaching a 
surface element i of a surface model of human body shape 
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is calculated according to formula (1) by inputting an 
amount of direct solar radiation Idn to the normal plane 
and a transmittance tj of solar radiation of each of wall 
members (translucent members) through which solar 
radiation passes. 



m : a case that the sun can be seen 

1 drill* j co% &i through an m-number of wall members 
1 di = \ j=l through which solar radiation passes 

,0: Other cases 

Formula (1) 



The transmittance of solar radiation and the 
10 reflectivity of a translucent member have incident angle 
characteristics as shown in Fig. 12. Namely, when the 
incident angle to the translucent member becomes large, a 
reflection component will increase and the transmittance 
of solar radiation of the translucent member will 
15 decrease. There causes error in the amount of the solar 
radiation heat unless the incident angle characteristics 
are not considered. 

2) Calculation of an amount of sky-diffused solar 
radiation 

20 An amount of sky-diffusion solar radiation I si 

reaching a surface element i of a surface model of human 
body shape is calculated by an amount of sky-diffused 
solar radiation I sky to the horizontal plane and an angle 
factor of a wall surface (translucent member) through 

25 which solar radiation passes according to formula (2) . 

1 si = (Z F ij t j^V ' 0 ' 91 ^sky Formula (2) 

J 
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where Fij is a angle factor between surface elements i 
and j and pjj is a flag determining whether or not the 
reached solar radiation is sky-diffused solar radiation. 
A numerical value of 0.91 is a coefficient based on the 
5 incident angle characteristics of the transmittance of 
solar radiation of a translucent member (a glass sheet) 
with respect to sky-diffused solar radiation. 

3) Calculation of an amount of solar radiation 
reflected on the ground 

10. An amount of solar radiation reflected on the ground 

undergoes influence of the shape, the reflectivity, the 
direction characteristics of the ground, and accordingly, 
it is difficult to obtain a correct calculation value. 
In this embodiment, in the same manner as the method for 

15 calculating the amount of sky-diffused solar radiation, 
an amount of solar radiation I gi reflected on the ground 
to reach each surface element of the surface model of 
human body shape is calculated by using an amount of 
total sky-diffused solar radiation I ho i to the horizontal 

20 plane and the albedo (the reflectivity of solar radiation 
energy) p g of the ground according to formula (3) . 

{ ,J gi =<Z F iJ t jrif'0.91) / ? g I ho2 
j Formula (3) 

Jhol = 1 dn sinh+ I sky 

where yij is a flag for determining whether or not the 
reached solar radiation is a light reflected on the 
25 ground, and h is a solar altitude. 



4) Calculation of an amount of relative reflection 
When solar radiation reaches a surface element of 

the surface model of human body shape, relative 
reflection is generated in response to the reflectivity 
of the surface of the surface model. The reflection 
includes a diffused reflection, a specular reflection and 
a complex reflection comprising both types of 
reflections. In this embodiment, a completely diffused 
reflection is assumed in order to simplify calculation. 
Under this condition, an amount of relative reflection 
can be calculated by a radiosity method using an angle 
factor between surface elements according to the formula 
(4) . 

j Formula ( 4 ) 

5) Method for calculating an amount of absorptance 
of solar radiation at the surface of a thermal manikin 

An amount of direct solar radiation, an amount of 
sky-diffused solar radiation, an amount of solar 
radiation reflected on the ground and an amount of 
relative diffused reflection calculated for each surface 
element of the surface model are all converted to an 
amount of absorptance of solar radiation on the surface 
of the surface model. The method of conversion to the 
amount of absorptance of solar radiation at the surface 
of the surface model is made by multiplying the above- 
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mentioned amounts of solar radiation by an absorptance of 
solar radiation ai of a surface element i of the surface 
model according to formula (5). 

I i = a i (I di + I si + I gi + I ri ) Formula (5) 

5 By treating this as an amount of absorptance of 

solar radiation (an amount of heat generated at the 
surface of a thermal manikin) in the equation of heat 
balance at a wall surface as shown in Table 7, it becomes 
possible to conduct a combined analysis of temperature 

10 and airflow. At the completion of the combined analysis 
of temperature and airflow, it is balanced with an amount 
of heat generated at the surface of the thermal manikin, 
an amount of thermal conduction by convection in response 
to an ambient air temperature of the thermal manikin, an 

15 amount of thermal conduction by radiation in response to 
the emissivity of the surface of the thermal manikin and 
an amount of total thermal conduction in clothes. 

Then, indoor physical quantities such as 
temperature, velocity and turbulence in the indoor space 

20 between the thermal manikin and the wall are calculated 
by utilizing a technique of CFD (Computational Fluid 
Dynamics) (Step 3 06) . Namely, various boundary 
conditions at surfaces of the thermal manikin and walls 
are determined, and a numerical simulation of the 

25 convection (including natural convection and forced 

convection) in the indoor space area is conducted for 
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each of the cubic elements to obtain the velocity, 
pressure, temperature, humidity and so on. 

Boundary conditions in terms of heat are previously 
stored in the external memory device 21 shown in Fig. 1. 
5 The external memory device 21 stores also thermophysical 
properties of solar radiation (the absorptance of solar 
radiation and the transmittance of solar radiation) at 
surfaces of interior part and exterior part, 
thermophysical properties on heat radiation (such as the 

10 emissivity) at surfaces of interior parts, a heat 

conductance between parts, the reference temperature at 
the exterior side and the convective heat transfer 
coefficient at the interior-side wall surface. 

As a concrete technique of CFD, a numerical analysis 

15 of a Navier-Stokes equation is made by using, for 

example, a finite element method, a finite volume method 
or a finite differential method. In particular, in this 
embodiment, a standard k-e model in a non-isothermal 
field is used. Then, the thermo-regulation reaction of 

20 the thermal manikin is calculated by utilizing the 
thermo-regulation model shown in Fig. 7 (Step 3 07) . 
Then, calculating operations (Steps 306 and 3 07) are 
repeated until the skin temperature (or the internal 
temperature) of each surface element of the surface model 

25 of human shape is converged to a predetermined value 
(Step 308). Then, the equivalent temperature is 
calculated at Step 309. 
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The combined calculating operations of CFD , the 
radiation and the thermo-regulation model can be obtained 
by solving the equation of heat balance as shown in Table 
7. Namely, a solution of the equation of heat balance 
5 comprising an amount of convective heat transfer at the 
wall surface of the passenger compartment, an amount of 
radiative heat transfer, an amount of absorptance of 
solar radiation and an amount of heat transfer at the 
exterior side is obtained, and calculations of 1) 

10 connection from the convection field to the radiation 

field, 2) convergence of the thermo-regulation model and 
3) connection from the radiation field to the convection 
field on the surface of clothing are conducted. In the 
calculation of 1) , the temperature Ti at the surface of 

15 clothing and an amount of radiative heat transfer Q r i<net) 
are obtained; in the calculation of 2), the calculation 
of the thermo-regulation model is conducted to obtain the 
temperature Ti at the surface of clothing and the skin 
temperature T re f are obtained, and in the calculation of 

20 3), the reference temperature T in at the indoor side and 
the temperature Ti at the surface of clothing are 
obtained. Further, the above-mentioned calculations may 
be conducted in consideration of humidity. 



Table 7 
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<Equation of heat balance at a wall surface> 

Amount of Amount of \ Amount of 
convective radiative J heat 
heat heat / transfer at 

transfer transfer I the outdoor 

\ side 

Amount of 

absorption of 

solar radiation 

/[— ] : Number of surface element 

T in [°C] : Reference temperature at the indoor side 
flfci[W/m^ *°C] : Convective heat transfer coefficient 

K '[W/m^ • °C] : Heat transfer coefficient between indoor 
side wall surface and the reference 
temperature at the outdoor side 

7,[W/m ] : Amount of absorption of solar radiation 
T re y[°C] : Reference temperature at the outdoor side 
Tj[°C] : Temperature of wall surface 
q ( we /)[W/m^] : Amount of radiative heat transfer 

( "Net " is a balance between a received radiation and 
a generating radiation.) 



<Heat balance at the surface of clothing> 

1) Convection field -> radiation field 

*ki(Tin -Ti) + Qr t {net) + h + K \T ref - T t ) = 0 

K '[W/m^ -°C] : Heat transfer coefficient of clothing 
Tj n [°C] : Reference temperature at the indoor side 
T re y[°C] : Skin temperature 

2) Convergence calculation of thermo-regulation model 

Same as 1) 

3) Radiation field Convection field 

Same as 1) 
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Particulars of calculating the thermo-regulation 
model in 2) are shown in Table 8. Here, the equation of 
heat balance is applied to all surface elements of the 
surface model of human shape. 



Table 8 



10 



15 



20 



25 



Skin layer 



Heat 
source 



'sk 



PC 



T sk =36.4-0.54.05* 



*1: Temperature sensor 



The equation of heat balance comprises Q sk as a 
sensible heat loss, 3 6.4°C as the inner temperature of 
thermal manikin and 1/0.054 (W/m 2 k) as the heat transfer 
coefficient of a portion from the position at which the 
inner temperature is determined to the skin surface. 
Accordingly, the calculation of the thermo-regulation of 
the thermal manikin is conducted by keeping a balance 
between the skin temperature T sk of the thermal manikin 
and the sensible heat loss Q sk . In this model, it is 
unnecessary to conduct the simulation of heat transport 
by perspiration or blood flow as required in the 65MN 
model, and the construction of a mathematical model can 
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be simplified to thereby reduce the amount of the 
calculation. Further, in order to simplify further the 
thermo-regulation of the thermal manikin, an amount of 
heat generation by a heat source may be adjusted so that 
the temperature of the skin layer of the thermal manikin 
is constant or an amount of heat generation by the heat 
source may be constant, or perspiration or breath may be 
considered. 

After the above-mentioned series of calculations are 
made in the operation part 1.1a in Fig. 1, the calculation 
of the equivalent temperature (Step 3 09) is made. Fig. 
13 shows a flow chart showing particulars of Step 309. 
First, 1) an amount of heat loss from the skin surface 
and 2) an averaged skin temperature are calculated from 
physical quantities of each part of human body (Step 
401) . Finally, the equivalent temperature at the skin 
surface of the thermal manikin is calculated by solving 
the equation shown in Step 402. 

In the flow chart of Fig. 13, the equivalent 
temperature of the whole body is calculated. However, 
calculation may be conducted for each surface element or 
each part of the surface model of human shape. In this 
case, an amount of heat loss and a skin temperature (or 
an inner temperature) are obtained for each surface 
element or each part of the surface model. Thus, by 
calculating the skin temperature and so on for each 
surface element, the temperature distribution in each 
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part can be simulated in detail and a correct thermal 
sensation can be calculated. These results of simulation 
are displayed in the display device 2 3 after these 
results are fabricated into a graph, a Table, a still 
5 picture (or a skin temperature distribution or the like) , 
animation or the like by the control of the evaluation 
part lib in Fig. 1. The performer of the simulation can 
determine the thermal comfort of a translucent member 
based on these displays. 
10 EXAMPLE 

Now, the present invention will be described in 
further detail with reference to Examples. However, it 
should be understood that the present invention is by no 
means restricted to such specific Examples. 

15 As a space to be analyzed, a vehicle body of RV 

(recreational vehicle) type is assumed. A result of 
simulation in a case that a passenger is on the driver's 
seat, is shown. 

Fig. 14 is a diagram showing the shape of a vehicle 

20 50. The vehicle accommodates a thermal manikin 51 as a 
passenger. A simulation was conducted under calculation 
conditions that the passenger compartment of the vehicle 
was subjected to cooling and air-conditioning in 
midsummer, date was August 31, the location was Tokyo, 

25 the time was 16:00 and the vehicle was running in a 
direction of around west at 40 km/hr at an outer 
temperature of 33°C. The thermal manikin 51 was assumed 
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to wear a typical summer clothing. 

In order to reproduce an outdoor condition in an 
artificial climate control room, 36 infrared lamps were 
disposed in front of the vehicle body 50 so as to provide 
5 the same amount of solar radiation to the horizontal 
plane in an outdoor space. Then, simulation was 
conducted by irradiating light beams of these lamps to 
the vehicle body 50. 

Table 9 shows a layout of glass sheets. Study was 

10 made on two assumed cases that green grass was used for 
all windows as standard equipment and a high-functional 
glass having a solar radiation controlling function was 
used for all windows. Table 9 describes kinds of colored 
glass, and thicknesses. In Table 9, data of each of two 

15 glass sheets constituting a lamination glass for the 
windshield are described. 



Table 9 





Windshield 


Front 
door 
glass 


Rear door 
glass 


Rear 
quarter , 
rear 
glass 


Case 1 
(Normal 
type) 


Green (2 mm) + 
green (2 mm) 


Green 
(3.5 mm) 


Green 
(3.5 mm) 


Green 
(3.1 m) 


Case 2 
(Solar 
radiation 
controlling 
type) 


Green (2 mm) + 
green (2 mm 
with coating) 


UV cut 
green 
(3.5 mm) 


Privacy 

gray 
(3.5 mm) 


Privacy 

gray 
(3.1 m) 



20 
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Fig. 15 shows an amount of solar radiation absorbed 
in the vehicle body and the surface of the human body in 
a case of using the solar-reduction glass. When a 
numerical simulation was carried out under outdoor 
5 conditions shown in Fig. 15(a), it was found that a chest 
portion of the thermal manikin absorbed solar radiation 
in an amount of as large as about 22 0 w/m 2 . On the other 
hand, when a numerical simulation was carried out under 
the condition of the artificial climate control room 

10 shown in Fig. 15(b), it was found that the chest portion 
of the thermal manikin absorbed a largest amount of solar 
radiation of about 70 w/m 2 . It is considered that this 
difference is caused due to a difference of light 
transmit tance even though the same kind of glass is used, 

15 because there are differences of spectral properties of 
light sources and spectral transmission characteristics 
of glass. 

Fig . 16 is a graph showing amounts of solar 
radiation absorbed in a human body. Under an outdoor 

20 condition, there is a difference of 15% in the amount of 
solar radiation absorbed in a human body due to using 
different types of glass. However, under the condition 
of the artificial climate control room, the difference is 
large as about 65%. 

25 Fig. 17 is a graph showing heat loads by air- 

conditioning. Under an outdoor condition, there is found 
a difference of about 5% in the heat load by air- 
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conditioning due to using different types of glass. 
However, under the condition of artificial climate 
control room, the difference is large as about 12%. 

Fig. 18 is a graph showing equivalent temperatures 
5 of typical 16 parts and the whole body of a passenger on 
the driver's seat. Equivalent temperatures within the 
acceptance limit means generally comfortable. It is 
understood that there is a difference of about 10-15°C 
depending on parts which receive solar radiation, and 

10 there is a large influence to thermal sensation because 
glass interrupt partially solar radiation. It is 
understood that the equivalent temperatures deviate from 
the acceptance limit under the outdoor condition even 
though either the glass case 1 or the glass of case 2 is 

15 used. Further, there is found a difference of about 1- 
3°C in view of the equivalent temperature of the whole 
body. It is understood that when the high-functional 
glass is used under the indoor condition, the equivalent 
temperature is slightly lower than the acceptance limit, 

20 and it provides a thermal sensation of "slightly cool". 

As is clear from the above-mentioned explanation, it 
was confirmed that the same result as the case using the 
conventional thermal manikin could be obtained by 
conducting a numerical simulation of the thermal model of 

2 5 human body 3 0 shown in Fig. 5 by using a computer, 
without conducting measurements using the thermal 
manikin . 
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Although explanation has been made as to an example 
using the glass sheets as a translucent member, the 
present invention is not limited to such example but is 
applicable to, for instance, an organic resin plate (a 
5 polycarbonate plate, an acrylic plate etc.) or an organic 
resin film usable for organic glass, vinyl houses and so 
on. Further, as the standard effective temperature for a 
human body, a standard effective temperature obtained by 
evaluating and calculating influence of each part more 

10 precisely than 6 5MN SET* may be used. 

Further, description has been made as to the 
simulation conducted in a steady state. However, it is 
understood that a simulation in an unsteady state can be 
carried out by applying the above-mentioned series of 

15 procedure for each predetermined time step. For example, 
when the thermal comfort of a glass sheet used for a 
vehicle is to be evaluated, it is considered that the 
time step is determined to be a range of from several sec 
to several min, and calculation is conducted over a range 

20 of from several min to several ten min (or several hr) . 
The translucent member in the present invention 
includes members to be attached to a structure to 
introduce light therein, such as a single glazing glass, 
an insulated glazing glass, a lamination glass comprising 

25 an organic resin film (a polyvinyl butyral film or the 
like) and a plurality of glass sheets sandwiching the 
organic resin film, an organic resin film or an organic 
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resin plate. The translucent member is not necessarily 
completely transparent. For example, glass colored by 
adding iron, cobalt or the like to absorb heat or a thin 
metallic film may be coated thereon to reflect heat wave. 
5 Further, as the translucent member, a flat plate or 

curved plate is generally used. However, it may have 
another shape as long as it has function to introduce 
light. Further, description has been made on the 
assumption that the humidity in the indoor space is 

10 constant in order to simplify explanation. However, it 
is clear that the humidity can be calculated for each of 
the cubic elements. Further, in order to simplify 
calculation, a simulation may be conducted under the 
condition that the temperature and airflow fields in the 

15 indoor space are constant. 

In addition, a simulation may be conducted in 
consideration of the outer configuration of the vehicle 
shown in Fig. 2 and ambient of the vehicle. By doing 
this, correct evaluation about influence of the outer 

20 configuration to the passenger's thermal comfort. 

As described above, in accordance with the present 
invention, the material property of each part capable of 
providing the optimum thermal comfort can be obtained by 
a numerical simulation of thermal comfort index of a 

2 5 measuring device having a shape imitating a human shape, 
located in an indoor space or a passenger compartment. 
In particular, when calculation of the thermo-regulation 
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of a thermal manikin is carried out by (a) adjusting an 
amount of heat generated from a heat source so that the 
temperature of a skin layer of the thermal manikin is 
constant, (b) making an amount of heat generated from the 
5 heat source constant, or (c) balancing the temperature of 
the skin layer of the thermal manikin with a sensible 
heat loss at the surface of the skin layer, the 
simulation of the thermo-regulation function can be 
conducted more easier than the case of using a high-grade 

10 thermo-regulation model such as the 65MN model or the 

like. Namely, it is unnecessary to conduct a simulation 
of the heat transport by perspiration or blood flow as 
required in the 65MN model, and the construction of a 
mathematical model can be simplified to thereby reduce 

15 the amount of calculation. However, the present 

invention does not exclude a high-grade thermo-regulation 
model such as the 65MN model, and such thermo-regulation 
model may be selected optionally depending on the 
performance of computer used or accuracy of a required 

20 simulation . 

Further, when the skin temperature is calculated for 
each surface element of a thermal manikin to calculate 
the thermal comfort index for each part of it, the 
optimum design of a structure can be made even though the 

25 indoor space or the passenger compartment is in an 
ununiform environment in terms of heat. 

Further, the thermal comfort of a translucent member 
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can be evaluated by only a numerical simulation, and 
accordingly, a labor of preparing a prototype of a 
structure or a vehicle can be saved. Namely, evaluation 
of various kinds of window glass sheets is possible by 
5 only simple operations of changing parameters of the 

absorptance of solar radiation, the emissivity of solar 
radiation and so on. Further, since it is unnecessary to 
prepare the prototype, the thermal comfort of the 
translucent member can be evaluated in a shorter time and 
10 at a low cost in comparison with the conventional 
technique . 

Further, the present invention can be employed in 
order to develop new material for a translucent member 
capable of providing the most excellent thermal comfort 

15 depending on the shape of an indoor space or a passenger 
compartment. Further, in the present invention, since 
the thermal comfort of a translucent member can be 
expressed by an objective index, i.e., a thermal comfort 
index, it is possible to explain comprehensively the 

20 thermal comfort of a translucent member to building 

manufacturers, vehicle manufacturers, air-conditioning 
system manufacturers, material manufacturers, and general 
users. Further, it is possible to provide programs of 
the above-mentioned simulation or to execute the 

25 simulation through internet. Accordingly, the present 

invention can provide new services for evaluation of the 
thermal comfort of a structure or a support for designing 
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a structure. 

The entire disclosure of Japanese Patent Application 
No. 2002-279695 filed on September 25, 2002 including 
specification, claims, drawings and summary is 
5 incorporated herein by reference in its entirety. 



